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Abstract 

We present the first measurements of the double spin asymmetries Ann and Ass at \fs = 200 
GeV, obtained by the pp2pp experiment using polarized proton beams at the Relativistic Heavy Ion 
Collider (RHIC). The data were collected in the four momentum transfer t range 0.01 < \t\ < 0.03 
(GeV/c) 2 . The measured asymmetries, which are consistent with zero, allow us to estimate upper 
limits on the double helicity-flip amplitudes 4>2 and (p4 at small \t\ as well as on the difference Aar 
between the total cross sections for transversely polarized protons with antiparallel or parallel spin 
orientations. 

PACS numbers: 13.85.Dz and 13.88. +c 
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I. INTRODUCTION 



In this paper we present the first measurements of the double spin asymmetries Ann and 
Ass m elastic scattering of polarized protons at RHIC at */s = 200 GeV and low t-range 
0.01 < \t\ < 0.03 (GeV/c) 2 . The experiment was carried out by the pp2pp Collaboration 



QflflQ. 



The 



the energies reached in previous 



at the RHIC-Spin complex as a part of the pp2pp experimental program 
oresent results complement our previously published results on the analysing power An 
4j. The center of mass energy exceeds by a factor of 
measurements of the elastic scattering spin parameters 

The kinematic region of our experiment is suitable for an investigation of the spin depen- 
dence of diffractive scattering — the dominant mechanism at high energies. In particular, 
the physical motivation of the present experiment is connected with two long-standing ques- 
tions. First, is s-channel helicity conserved in elastic scattering at asymptotic energies, 
which is directly related to the question if the Pomeron exchange contributes to the spin 
effects (ref. jl0| and references therein)? Second, does the hypothetical Odderon exist 111]? 

The double spin asymmetry A NN is defined as the cross section asymmetry, 
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for both beams fully polarized along the unit vector h normal to the scattering plane. 
The asymmetry Ass is defined analogously, but for both beams fully polarized along the 
unit vector s in the scattering plane and normal to the beam. Elastic pp-scattering is 
described by five independent helicity amplitudes: two helicity conserving ones <p\ and 6^, 



io| 



two double helicity-flip ones 4> 2 and 04, and one single helicity-flip amplitude 5 (see Ref. 
for definitions). Each amplitude consists of hadronic and electromagnetic parts; <pi = 0f m + 
$l ad . The double spin asymmetries and Aoy = cr^ — <r^, which are the subject of this 
paper, relate to the helicity amplitudes as follows [T" 
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In the small t region the parameters Ann and Ass are sensitive to the interference 
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between hadronic spin-flip amplitudes cj>2 ad and <p ! l ad , and the electromagnetic non-flip am- 
plitude. This provides a sensitive tool to study the spin dependence of diffractive scattering 
at asymptotic energies and to search of the hypothetical Odderon exchange [llj|. Because 
Pomeron and Odderon have opposite C-parities, it is expected in leading order, that if 
Pomeron and Odderon have the same asymptotic behaviour, up to logarithms, they are out 
of phase by approximately 90° [12]. Therefore, if they couple to spin, their interference with 
the electromagnetic non-flip amplitude will result in different t dependences of the double 



spin asymmetries (for a more detailed discussion of the A NN case see [111]). 



II. THE EXPERIMENT 

The experiment pp2pp is located at the "2 o'clock" position of the RHIC ring. The 
layout of the experiment is shown in Fig. [TJ The two protons collide at the interaction 
point (IP) and since the scattering angles are small, the scattered protons stay within the 
beam pipe until they reach the detectors. The measured coordinates are related to the 
scattering angles by the beam transport matrix. The coordinates are measured by silicon 
microstrip detectors (SSD) positioned just above and below the beam orbits by insertion 
devices - Roman Pots (RP) [llj. Each RP contains four planes of SSDs (two vertical and 
two horizontal) to provide redundancy for the track reconstruction. High quality SSDs 
manufactured by Hamamatsu Photonics with sensitive area of 4.5 x 7.5 cm 2 and pitch of 
100 /im were used. The identification of elastic events is based on the collinearity criterion, 
hence it requires the simultaneous detection of the scattered protons in a pair of RP detectors 
on either side of the IP. More details on the experiment and the technique used can be found 
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In addition, the inelastic event rates were monitored by eight scintillation detectors, four 
on each side of the IP, positioned outside the beam pipe (cf. Fig. 1). To reduce background 
due to the beam halo the outer scintillators, labelled IP3, IP4 for the detectors on the side 
of the IP towards RP1, and IN3, IN4 on the opposite side, were used for this analysis. These 
detectors covered a pseudorapidity range of 4.2 < rj < 5.3. The inelastic trigger was defined 
as (IP3 OR IP4) AND (IN3 OR IN4) in coincidence with the beam crossing signal from 
RHIC. 

In order to prevent inelastic events from dominating the event rate and causing excessive 
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FIG. 1: Layout of the PP2PP experiment. Note the detector pairs RP1, RP2 and RP3, RP4 lie 
in different RHIC rings. Scattering is detected in either one of two arms: Arm A is formed from 
RP3U and RP1D. Conversely, Arm B is formed from RP3D and RP1U. The coordinate system is 
also shown. 

dead time, the inelastic triggers were prescaled by a factor of 500. The data were recorded 
in an event-by-event mode, and broad offline cuts were made on the TDC and ADC to select 
events, yielding approximately 1400 inelastic events for each of 50 bunch-pair crossings. 

Each bunch in one RHIC ring (B) collides with a bunch in the other ring (Y) to form dif- 
ferent spin combinations. In this experiment, 13 different bunch-pairs yielded collisions with 
B^Y^ bunches, 13 yielded B^Y-^ bunches, 11 yielded B^Y^, and 13 yielded B^Y^ combinations 

The inelastic event rate was used as one of the two methods to monitor the relative 
luminosity for collisions of bunches with different polarizations. The product of intensities 
of colliding bunches was used as a second method to monitor the relative luminosity. Both 
methods are described and compared in Chapter IV. 

III. SELECTION OF ELASTIC EVENTS 

Details on the identification of hits in the silicon detectors, hit clustering and selection of 
elastic events can be found in Ref. [4j. Because of the collinearity of the scattered protons 
a correlation between coordinates measured on each side of the IP is required. Hence, the 
main criterion to select the elastic scattering events was the hit coordinate correlation in the 
corresponding silicon detectors on the opposite sides of the IP. An example of the distribution 
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FIG. 2: Distribution of the sum of x-coordinates in a pair of silicon detectors on opposite sides of 
the IP. The coordinates are given in detector strip numbers. 

of the sum of x-coordinates measured by a pair of detectors on both sides of the IP is shown 
in Fig. [2j One can observe that the background under the elastic peak is small. To select 
an elastic event, a match of hit coordinates (x,y) from detectors on the opposite sides of the 
IP was required to be within 3a for x- and y-coordinates. 

The average detector efficiency was 0.98, and the upper bound of the elastic events loss 
due to all criteria was < 3.5%. 

The background originates from particles of inelastic interactions, beam halo particles and 
products of beam-gas interactions. The estimated background fraction varies from 0.5% to 
9% depending on the ^-coordinate. Since in our analysis the coordinate area was essentially 
limited to y > 30 strips, the background in the final sample does not exceed 2%. 




IV. DETERMINATION OF THE DOUBLE SPIN ASYMMETRIES 

After the above cuts, a sample of 2.3 million elastic events was collected in the interval 
0.010 < — t < 0.030 (GeV/c) 2 . The azimuthal angle dependence of the cross section for the 
elastic collision of the vertically polarized protons is given by 

27r— °- = • (1 + (P B + P Y )A N cos0 + P B P Y {A NN cos 2 + A ss sin 2 <p)) , (3) 
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where Pb and Py are the beam polarizations, Ann-, Ass ar e the double spin asymmetries 
and An is the single spin asymmetry. The double spin raw asymmetry 5(0) is 

5(0) = P B Py (A N n cos 2 + A ss sin 2 0) 

jytt (0) /£TT _|_ #U (0) /LH _ jy-ti (0) /L u - JV^ (0)/L iT 
~ ATTT(0)/LTT + ATH(0)/LH + ATU(0)/LU + JViT(0)/LlT ' ( > 

where L* J is the relative luminosity for the sum of bunches with a given spin combination 
{hj} £ (t1\ Tl) IT) ID- The raw asymmetry 5 was calculated as a function of the azimuthal 
angle using 5°-bins in the three t-intervals 0.010 < —t < 0.015, 0.015 < — t < 0.020, 
0.020 < -t < 0.030, as well as in the combined interval 0.010 < -t < 0.030 (GeV/c) 2 . 

To evaluate 5(0) two estimates of relative luminosities L lJ were used. For the first one 
the rates of 'inelastic' counters 

aa 

, Nlf iel , summed over bunch crossings with a given spin 
combination were used, L l £ el ~ Yli-^lnei- ^ s ^ ne secon d estimate the sums of the machine 
bunch intensities products for each spin combination L l ^ nt ~ ^ ' ly were used. Because 
the relative statistical errors and relative errors on the luminosity are comparable, they have 
to be combined for the determination of errors of raw asymmetries. 

An example of the double spin raw asymmetry 5(0) distribution for the whole t— interval 
is shown in Fig. [3] together with the fitted function P B Py (A N n cos 2 + A S s sin 2 0). A N n 
and Ass are the fit parameters and Pb ■ Py = 0.198±0.064. The displayed raw asymmetry 
was obtained using bunch intensities for an estimate of the relative luminosities. 

In order to facilitate separation of contributions of the helicity amplitudes 02 and 04 
to the double spin asymmetries (cf. Eqs 2), we performed also alternative fits to 5(0) = 
Pb Py (ai + a 2 cos 20) using a\ = (A N n + Ass)/2 and a 2 = (A N n — ^ss)/2 as fit parameters. 

As a cross check the 5(0) distributions were obtained using each of the two methods to 
estimate the relative luminosities and the results of the fits to the 5(0) distributions were 
compared. Both methods gave consistent results, although with larger errors for normaliza- 
tion by the numbers of inelastic events. Thus from now on we only present results obtained 
by using the bunch intensities. 

V. SYSTEMATIC ERRORS 

A detailed study of the errors on the relative luminosity estimated using bunch intensities 
was carried out. In the first stage of data processing the stability of the normalization was 
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FIG. 3: The raw double spin asymmetry 5{cj)). 

checked by inspecting for each spin state {i,j} the ratio N l ^ ast /Il ■ V Y as a function of the 
bunch crossing number. The bunch crossings for which a significant deviation of the ratio 



from the average was observed were excluded from the analysis. Next the L 



hJ 
bint 



^2 Ib'Iy summed over runs were evaluated for each position of detectors. The ratios between 
any pair of P J 's are quite stable as a function of the position of detectors, at the level of 
< 0.2%. 

For the method using bunch intensities to be applicable, it is mandatory that L\l nt = 
k ■ P' J = k ■ Ib ' ly an d the coefficient k not to depend on the spin state. The stability 
of Nlf ast / P'i was analysed on a bunch-by-bunch basis. The main result of this analysis is 
that the relative error on the normalization is equal to a(k)/k = 0.0124 ± 0.0015. From the 
analysis of the spread of k on a run-by-run basis we conclude that about half of the error 
on k is due to the accuracy of the intensity measurements. 

The systematic errors on the double spin asymmetries depend on the errors of the posi- 
tions of the beams in the x and y planes and on the transport matrix elements, especially 
the effective lengths. The systematic errors due to these two effects were estimated using 
the experimental sample and assuming Ax = 0.5 mm, Ay = 1.0 mm, AL X ^/L X ^ = 5% and 

ALV eff/ LV eff = 5% • 

The effect of these two sources on the systematic error depends on a particular asymmetry; 
the uncertainty on the effective lengths is the dominant contribution to the systematic error 
of Ann while the uncertainty on the beam coordinates is dominating the one for Ass- The 
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combined systematic errors for various double spin asymmetries are given in Table I. The 
systematic errors are about half of the statistcal errors (including the relative luminosity 
errors) . 

The polarization values of the proton beams were obtained from the accelerator po- 
larimetry [14J]. They were evaluated using An measurements for elastic proton-carbon (pC) 
scattering at small |£|-values, in the range 0.01 — 0.02 (GeV/c) 2 . The details are described 
in Ref. [9j. During the period in 2003 when the present data were taken the beam po- 
larizations were Py = 0.411 ± 0.021 and Pb = 0.481 ± 0.027. These values correspond 
to the present analysis, which includes bunches with all spin combinations. The indicated 
errors are purely statistical. In addition there is a systematic uncertainty of 13%, due to 
the calibration of the pC polarimeter, which is correlated for both beams. This gives for the 
combined statistical and systematic errors of the measurement of the product of the polar- 
izations Py ■ Pb — 0.198 ± 0.064. Thus, due to the uncertainties on the beam polarization 
the scale error on the measured asymmetries is equal to 32% . 



VI. RESULTS AND CONCLUSIONS 



The results on the double spin asymmetries for the whole t-interval 0.010 < — t < 0.030 
(GeV/c) 2 , at an average < —t >= 0.0185 (GeV/c) 2 , are presented in Table I. The most 
accurately determined asymmetry is Ass — 0.0035 ± 0.0081, which is consistent with zero 
at la confidence level. The asymmetry Ann = 0.0298 ± 0.0166 as well as the combinations 
{A NN + A ss ) /2 = 0.0167 ±0.0091 and (A NN - A ss )/2 = 0.0131 ± 0.0096 are also small and 
consistent with zero. 



TABLE I: Double spin asymmetries Ann, Ass, {Ann + Ass)/2 and {Ann — Ass)/2 for the 
t-interval 0.010 < -t < 0.030 (GeV/c) 2 at < -t >= 0.0185 (GeV/c) 2 . 





Ann 


Ass 


{Ann+A ss )/2 


{Ann - A ss )/2 


Asym 


0.0298 


0.0035 


0.0167 


0.0131 


AAsym (stat.+norm.) 


±0.0166 


±0.0081 


±0.0091 


±0.0096 


A Asym (syst.) 


±0.0045 


±0.0031 


±0.0034 


±0.0072 


AAsym due to A{P Y ■ P B ) 


±32.3 % 
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At collider energies one expects [lCj the two helicity conserving amplitudes 0i and 03 
to be equal, 0i ~ 03 ~ 0+ = (0i + 03)/2. In addition, if for simplicity one assumes the 
Coulomb amplitude to be pure real, 0+ m = — Im 0+ -t c /t, where at our energy t c = —0.0013 
(GeV/c) 2 , the general formulae for the transverse double spin asymmetries (cf. Eq. 2) 
reduce to approximate ones: 

A NN ■ lm 2 + = 2|0 5 | 2 + Im0 + {lm0 2 - Im0 4 + (p - t c /t) • (Re0 2 - Re0 4 )} , 

A S s ■ lm 2 + = Im0 + {lm0 2 + Im0 4 + (p - t c /t) ■ (Re0 2 + Re0 4 )} , 

{A NN + A ss )/2-lm 2 (f) + = |0 5 | 2 + Im0 + {rm02 + (p-t c /t)-Re0 2 }, 

(A NN - A ss ) 1 2 • lm 2 0+ = |0 5 | 2 - Im0+{lm0 4 + (p - tjt) ■ Re0 4 } , (5) 

where p = Re0+/Im0 + . For brevity, in Eqs[5]and in the following we omitted the super- 
scripts 'had' for hadronic amplitudes. 

As seen from the above formulae, taking (A N n + Ass)/2 and (A N n — A S s)/2 combinations 
allows one to select the amplitudes 02 and 4 , respectively. Let us first assume that the 
phases of 02 and 4 are about the same as the phase of the non-flip amplitude + . Then the 
last term in the brackets can be neglected and the la upper limits for the imaginary parts of 
02 and 4 at t = -0.185 (GeV/c) 2 are: Im0 2 /Im0 + < 0.028 and -Im0 4 /Im0+ < 0.026, or 



in terms of the conventional ratios r2 = 02/(2 Im0 + ) and r 4 = — m 2 4 /(tIm0 + ) lOj they 
are Imr 2 < 0.014 and Imr 4 < 1.25. 

A more precise limit on Im 02 at t close to zero and therefore on Acxr = crJot ~~ °tlt can be 
obtained using the t-dependence of the asymmetry Ass and extrapolating Ass to t = —0.01 
(GeV/c) 2 where the term containing the real parts of amplitudes vanishes. For that purpose 
the corresponding experimental distributions 5(0) in the three t-intervals were fitted with 
a function Pb Py {Ass + (^4 aw — Ass) cos 2 0} with Ass as a fit parameter. Here the term 
Ann — Ass was not fitted, but calculated as a predefined function of t. At small \t\ this 
term is proportional to t, A N n — Ass = C -t, because of the kinematical factors in 5 and 4 
resulting from angular momentum conservation n]] . The constant C was calculated using 
the value of Ann — Ass from the 5(0) fit for the combined t-interval (cf. Table I). The 
results of the fit of Ass in bins of t are given in Table II. 

With the linear extrapolation to t = -0.01 (GeV/c) 2 we obtain A ss (t ) = 0.0037 ± 
0.0104. Neglecting the contribution of 4 to Ass and the variation of 02 over the small 
range of t one obtains Im0 2 /Im0 + = 0.0037 ± 0.0104, Imr 2 = 0.0019 ± 0.0052 and Aa T = 
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TABLE II: Double spin asymmetry A$s m different intervals of t. 



-t interval (GeV/c) 2 


0.010 - 0.015 


0.015 - 0.020 


0.020 - 0.030 


— <t> (GeV/c) 2 interval 


0.0127 


0.0175 


0.0236 


A SS 


0.0005 


0.0076 


0.0015 


A A ss (stat.) 


±0.0071 


±0.0056 


±0.0061 



—0.19 ± 0.53 mb. The quoted errors are the quadratic sum of the statistical, normalization 
and systematic errors and the error due to the beam polarizations uncertainty. 

The existing experimental data on Acx^ at low energies are fairly well described by a 

n 

Regge model with cuts [lSj. Our result is consistent with the prediction of the model for 
Aa T at v/i = 200 GeV. 

As seen from Eq. 5 the asymmetry Ass depends on the sum of the amplitudes 02 and 
04 and this sum does not couple to the leading and subleading Regge poles (Pomeron, 
Odderon, p, u, /, a 2 ) [10]. Thus one can expect that at yfs = 200 GeV only Regge cuts 
may contribute. We consider the effect on Ass of a possible contribution of the Pomeron- 
Odderon cut exchange in the t-channel as discussed in Refs 11] and jit] ]. In case of such 
exchange the phase of the 02 amplitude is shifted by 90° relative to the amplitude 0+, and 
Im0 2 = — pRe02 and thus Ass ~ — t c /t ■ Re02/Im0 + . Using the Ass value at t — —0.185 
(GeV/c) 2 (Table I) one obtains Re0 2 /Im0 + = -0.050 ± 0.130 or Rer 2 = -0.025 ± 0.065. 
Though this value is well consistent with zero it leaves wide room for a possible Pomeron- 
Odderon cut contribution. 

Theoretical predictions for double-spin asymmetries in elastic proton-proton scattering 
at high energies and small momentum transfers have been recently presented in Ref. 16]. 
The magnitudes of Ann and Ass have been estimated using results from an earlier deter- 
mination of the spin-couplings of the leading Regge poles [17j and the required Regge cuts 
were estimated using the absorptive Regge model. As the Odderon spin coupling is totally 
unknown, the predictions are given for various assumptions: (a) no Odderon, (b) weak Odd- 
eron spin coupling - equal to that of the Pomeron, (c) strong Odderon spin coupling - equal 
to the p Reggeon spin coupling. For none or a weak Odderon coupling the predicted values 
of the A NN and A S s asymmetries are very small. At y/s = 200 GeV and 0.01 < \t\ < 0.03 
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(GeV/c) 2 their values are in the range 0.001 - 0.002. On the contrary, for a strong Odderon 
spin coupling (like p) the double-spin asymmetries become significantly larger, at least by a 
factor of 10. Our results on the t-dependence of A S s (cf. Table II) support predictions of 
Ref. [l6] which assume none or a weak spin coupling of the Odderon. 

In conclusion, these are the first measurements of the transverse double spin asymmetries 
and the first results on the double helicity-flip amplitudes in the small \t\ region in elastic pp 
scattering at collider energies. From the measured double spin asymmetries we determined 
the parameters Imr 2 = 0.0019 ± 0.0053 and Aa T = —0.19 ± 0.53 mb, both being consistent 
with zero within errors. We also estimated the upper limit on Imr4 which is Imr4 < 1.25. 

Assuming the Pomeron- Odderon cut exchange one finds Rer 2 = —0.025 ± 0.065. The 
signs and central values of the real and imaginary parts of r 2 agree with expectations for 
Pomeron- Odderon cut exchange. Their magnitudes are consistent with an assumption of 
about 5% ratio of the cut amplitude to the dominant one. 
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